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Abstract. The Ising exchange interaction is a limiting case of strong exchange 
anisotropy and represents a key property of many magnetic materials. Here we find 
necessary and sufficient conditions to achieve Ising exchange interaction for metal sites 
with unquenched orbital moments. Contrary to current views, the rules established 
here narrow much the range of lanthanide and actinide ions which can exhibit Ising 
exchange interaction. It is shown that the Ising interaction can be of two types: (i) 
coaxial, with magnetic moments directed along the anisotropy axes on the metal sites 
and (ii) non-coaxial, with arbitrary orientation of one of magnetic moments. These 
findings will contribute to purposeful design of lanthanide and actinide based materials. 
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1. Introduction 

Strong magnetic anisotropy on the metal sites gives rise to novel magnetic properties, 
such as single-molecule magnet behavior ED 0. magnetic multipole ordering m m, 
and various exotic electronic phases [51 ©• In lanthanides and actinides the spin- 
orbit coupling exceeds the crystal-held splitting of the ground ionic LS term leading to 
unquenched orbital momentum L [7|. As a result the low-lying spectrum on these metal 
ions is well described as crystal-held split eigenstates of the total angular momentum 
J = L + S, where S is the spin of the corresponding term. If the symmetry of metallic 
sites is lower than cubic, the ground crystal-held multiplet is either a Kramers doublet 
or a quasi doublet, for odd and even number of electrons on the metal ion, respectively. 
The presence of unquenched orbital momentum in these doublets makes them strongly 
anisotropic, which is rehected also in the strong anisotropy of their exchange interaction. 
The limit for this anisotropy is the exchange interaction of Ising type, which in the case 
of two interacting doublets or interacting doublet and isotropic spin has the form: 

Hi J7” s i S'2z2 ■ 
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where s t is the pseudospin 1/2 describing the doublet state of site i, S 2 is the isotropic 
spin of site 2 and Zi denotes the projection of pseudospin or isotropic spin along the 
main magnetic axis of doublet on site i. 

The common beliefs concerning the Ising exchange coupling in lanthanides and 
actinides can be summarized in the following rules: 

1) . The interaction involving quasi doublets of metal ions with even number of 
electrons (the so-called Ising ions [3j) will always be of Ising type. However, in the 
case of metal ions with odd number of electrons the Ising interaction is only achieved 
when the corresponding Kramers doublets have zero magnetic moment in the transversal 
directions with respect to the main magnetic axis (i.e., are perfectly axial ®) 0- 

2) . I 11 both above cases the Ising interaction is of coaxial type (JTJ) , when the local 
magnetic moments in the exchange states are directed either along the main magnetic 
axes Zi of the corresponding metal sites or, in the case of isotropic spin, along the main 
magnetic axis of the neighbor site (Fig. 0a)) 0. 

Rule 1) appears to be weakly restrictive predicting, in particular, Ising exchange 
interaction to all metal ions with even number of electrons. As discussed below, both 
these rules are based on an oversimplified treatment of exchange interaction for metal 
ions with unquenched orbital momentum. Given the increasing interest for strongly 
anisotropic magnetic materials, the knowledge of precise conditions for the realization 
of Ising exchange interaction between strongly anisotropic metal ions would be highly 
desirable. 

In this work, the general conditions for achieving Ising exchange interaction for 
metal sites with unquenched orbital momentum are found via analysis of microscopic 
exchange Hamiltonian. These conditions are specified for interacting anisotropic metal 
ions in their doublet states and for such ions interacting with isotropic spins. Despite 
the overwhelming complexity of exchange interaction between centers with unquenched 
orbital moments, the realization of Ising exchange interaction is shown to depend in 
most cases solely on the structure of doublet wave functions on the metal sites and 
not on the details of intersite interaction. This reduces the task of designing magnetic 
materials with Ising exchange interaction to creation of appropriate crystal field on the 
metal sites. 

2. Exchange interaction for J-multiplets 

The ground-state (quasi) doublet wave functions in lanthanides and actinides can be 
written as linear combination of eigenstates of the corresponding atomic J-multiplet: 


j 


| M) = Y, C m \Jm), 


m=—J 


- m ) = e\ m), 


( 2 ) 
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where 6 is the operator of time inversion [lOj. The composition of these wave functions 
is sensitive to the details of the local crystal field and can involve, in particular, all 
eigenfunctions | Jm). This means that an adequate description of exchange interaction 
for these doublet states should involve an interacting Hamiltonian acting on the entire 
ground J-multiplet, i.e., expressed in terms of total angular momentum operator, J a , 
a = x,y,z. The latter was widely supposed to be of the Heisenberg-like form [3]: 

Htieis = JJl ‘ 32- (3) 

Despite the lack of justification, this form is often used for the description of interaction 
between lanthanides or actinides (or a similar form, oc Ji ■ S 2 , in the case of their 
interaction with an isotropic spin) [HI 121 Es El US, El El IE] - Direct calculations 
show that with this form of exchange interaction the rules 1-2) hold true [|| This is 
easily seen when we choose a particular form of wave functions ([2]) for strongly axial 
doublets: 

| ± Mi) = | Ji ± mi), i = 1,2, (4) 

achievable either at high axial symmetry [8] or at strong axial component of the crystal 
field on sites [2Uj. Indeed, it can be checked directly that for odd numbers of electrons 
the transverse components of any angular momentum is zero for m,; > 1/2, while in the 
case of even number of electrons the smallest rrii for a quasi doublet 0 is m j = 1 (the 
same for Eq. 0) . Under such circumstances the rules 1-2) and the conditions for zero 
transverse magnetization on metal ions (perfect axiality) hold true simultaneously [§] 
One should note that magnetic dipolar interaction between magnetic sites is also 
bilinear in total angular momentum operators albeit has a form different from Eq. ([3]): 

Hdip = - [Ji • J 2 - 3(J, ■ n 12 )(J 2 ■ n 12 )], (5) 

-^12 

where g^ is the Lande g factor, /i B is Bohr magneton, J? 12 is the distance between 
magnetic sites, and ni 2 is the normalized direction vector from site 1 to site 2. Therefore, 
it will reduce to Eq. (|TJ) too when downfolded on the axial doublets 0 K > 1 / 2 ). 
The corresponding constant J7di p in 0 is a simple expression of relative orientations of 
local main magnetic axes and the vector connecting magnetic sites |9j. This interaction 
is considered dominant for pairs of lanthanides ions, which justifies the description 
of their interactions via bilinear operators 0 [El [121 EE3 El 13 El El IE.- as 
well as the validity of rules 1-2). Such a view is, however, ungrounded for simple 
reasons that (i) dipolar magnetic coupling can be diminished arbitrarily on geometric 
grounds and (ii) the exchange interaction can significantly exceed the dipolar coupling, 
especially, when it is mediated by strongly covalent bridging groups [22]. In the cases 

f The same is true for the phenomenological Lines model m, H-Lines = USi ■ S 2 , where Si and S 2 are 
the spins of the ground state terms on the corresponding sites. 

§ A simplified treatment of kinetic exchange interaction within so-called 1/C7 approximation [21] gives 
for two interacting doublets |4| on Dy 3+ ions an Ising exchange interaction similar to the first Eq. (J0 . 
At the same time the interaction of one such Dy 3+ doublet with an isotropic spin results in an Ising 
exchange, as in the second equation of (JTj), only for mi = Ji = 15/2 (unpublished results). 
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Figure 1. Two types of arrangement of local magnetic moments by the Ising exchange 
interaction: type I - coaxial (a) and type II - non-coaxial (b) . The main magnetic 
axes on metal sites (dashed lines) are generally non-coplanar. Dotted lines indicate 
the main magnetic axis of the neighbor site. 


of lanthanides interacting with isotropic spins (transition metals or radicals) and of 
interacting actinides (between themselves or with isotropic spins), the magnetic dipolar 
interaction is significantly smaller than the exchange interaction. The latter, however, 
is not described by the simple form (|3| but represents a much more involved expression. 
In the case of two interacting ions with unquenched orbital momenta it has the form 

0123 : 


H = 


^ ' \Jkqk'q' 

kqk'q' 


QjPQohn) 

o»(Ji)o»,(j 2 y 


( 6 ) 


where Jkqk'q' is the exchange parameter, O^(Jx) and 0\,{ J 2 ) are Stevens operators [jj 
whose ranks k and k' have to obey the relation k + k' = even due to the invariance of 
the Hamiltonian with respect to time inversion m- In Eq. (|6]), the terms such that one 
of k and k' is zero are not included because they are not exchange but crystal field. A 
metal ion with unquenched orbital momentum interacts with an isotropic spin as follows 


H — ^ ^ yJkqlq' 


kqq' 


Og(Ji)S 2g - 

0»(J,)S 2 ’ 


(7) 


where k is odd due to the time reversal symmetry. 

The highest rank k of Stevens operators entering Eqs. 0 and (J7|) is found from 
the relation [23J: 


Lx = min [lli + l,2Ji], 


( 8 ) 


where l\ is the atomic orbital momentum and J\ is the total angular momentum of the 
ground atomic multiplet on site 1 (a similar relation holds for k' max on site 2). In the 


The Stevens operator is defined so that it transforms as spherical harmonics under rotation, which 
enables us to apply Wigner-Eckart theorem (see e.g., |Mj)- 
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case of f N ions, fc max = 21 + 1 = 7 for N = 2-4 and 7-13, fc max = 2J = 5 for N = 1, 5, 
and fc max = 0 for N = 6. 

The maximal component q (> 0) in Eqs. (J6]) and 0 is generally determined by 
details of intersite interaction [23]: 

Qmax min [A max + 1, /c max ] , (9) 

where A^ ax is the maximal difference of orbital momentum projections on site 1, m and 
n, entering as indices in the parameters V mrn ’ n ' n and t^ m ,t 2 ^ nl defining the direct and the 
kinetic intersite exchange interaction, respectively ( V is the bielectronic integral and t 12 , 
t 21 are electron transfer parameters). For low-symmetric systems, all these parameters 
will be non-zero, leading to A t ' nax = 2l x . Therefore, in such cases g max is determined by 
intrasite properties only (as /c max ). The maximal value of q' at the second site, q' max , is 
determined in a similar way. 

With the above information on the exchange interaction for entire J multiplet, the 
effective Hamiltonian is derived. For more information on the exchange interaction, 
such as the form of Jkqk'q 1 and explanations of & max and g max , see Ref. [23]. 


3. Derivation of the effective exchange Hamiltonian 

In this section, we derive the effective exchange Hamiltonians H e ff between the ground 
doublets of metal sites and between the ground doublet of the metal site and isotropic 
spin. This is done by projecting the complete form of the exchange Hamiltonian for the 
entire J-multiplet, Eqs. © and ((rj) , into the one within the truncated Hilbert space, 
and then describe the doublets by the s — 1/2 pseudospin operators. We consider first 
doublets of type Q, for which the derivation is the simplest. 


3.1. Interacting axial doublets 


Within the truncated Hilbert space, {\Jihi J 2 // 2 ) : hi = the exchange Hamiltonian 

0 reduces to 

H c ff = P12HP12, ( 10 ) 


where Pi 2 = P 1 P 2 , and P ?; = (* = 1,2) is the projector into the 

truncated space. The matrix element of the exchange Hamiltonian (Jo]) is 


EE Jkqk'q 1 


kk' qq' 


s~iJ 1 M 1 /^<J2k'2 

^ Jl/ilkq^ J2[12k' q' 
JlJlkO J2J2A/O 


(ii) 


as 


using the Wigner-Eckart theorem and C^ b p is the Clebsch-Gordan coefficients [25] ^j| 
Thus, substituting Eq. (JTTj) into Eq. (10), the effective Hamiltonian is obtained 

even ^ 

tt _ \ ^ \ ^ Ukqk'q' 

/ -j / -j J 2 J 2 

kk' qq' JiJifcO 


% We follow the convention of the spherical harmonics and the Clebsch-Gordan coefficients of Ref. [25] . 
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x [CjimifctA.o (| + |Ji - - mi|) 

+ (^,2m 1 |^imi)(J 1 - 77li| + <5 9 ,_2 mi | Jl “ mi)(Ji77l 1 |)] 

X (| ^ 2 ^ 2 ) (^2^21 + | ^2 “ m 2 ) (J 2 - m 2 |) 

+ (V,2ma I >^ 2 ^ 2 ) (</ 2 “ "U| + V ,-2m 2 |^2 ~ ^ 2 ) (^^l)] 


odd 


EE 


tJkqk' q' 


/^iJlJi ^ J 2-/2 
kk' qq' JiJikO J 2 J 2^'0 


x [C , S 1 1 fc0 ^ i o(|Jimi)(Jimi| - | Ji -?«!)(J x -mi|) 

+ - mil -^-arrnl^i -m 1 )(j 1 m 1 |)] 

x [C|r 2 \-'o^,o {\J 2 m 2 )(J 2 m 2 \ - | J -2 ~ m 2 )(J 2 -m 2 1) 

+ ( V,2m 2 | ^2'"U) (^2 “^ 2 ! ~ ^',-2m 2 |^2 ~ ^ 2 ) (^2'"U1)] ■ 


( 12 ) 


Here, the property of Clebsch-Gordan coefficient related to the time inversion is used 


Caab/3 ~ ( 1) 


a+b—c^yc —7 


a—a.b—( 3 ' 


(13) 


The two doublet functions (jij) on site z can be put in correspondence to the 
eigenfunctions 11/2, ±1/2) of the pseudospin 5* = 1/2: 


Ji^ij | 2^ "b 2^ $i z ii (14) 

Ji TTii)(Ji -/* ( 15 ) 

- mil = ~ v^Sj+i = s i+ = s iXi + is iyi , (16) 

Ji r ^i) (J%m i I V ^^2 —1 S^ — ^^iya ( 17 ) 


where i = 1,2, I % is the two dimensional unit matrix, s x ± = s lxi ± is\ yi and 
s 2 ± = s 2x . 2 ± is 2y2 ] Xi and y t denote Cartesian axes of site i, which generally do not 
coincide with similar axes of another site. Using the 1/2-pseudospin operators, 


Heff - X] 5 Z 


^Tkqk'q' 


^ r^JiJi rtJ2^2 
kk' qq' JiJifcO 


X 

X 


C jlmlk^q,oh + (^q, 2 mjl+ + dq,- 2 mi Sl-) 

Cjlmlk'^q'flh + C J2-m 2 k'2m2 (^',2m 2 S2+ + &q',-2m 2 S 2 -) 


odd 


EE 


<Jkqk' q' 


riJl Jl r^J2J2 
kk' qq' JiJikO J 2 J 2 WO 


X 

X 


(<5,,2 mi S.+ - Va™*i-) 

Si 

C&oVoi® + C&,«™ (V,2m,S 2+ - V,-2™ a S 2 _) 

S2 


(18) 
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Finally, the general exchange interaction (|6]) projected onto doublet functions of two 
sites can be recast in the following pseudospin exchange Hamiltonian: 


H e fr 71111\ l~2 4Kl zz Si Zl S2z2 

+ 7 Ti+i\ 1S2+ + Ji-hh- + 7J+1S1+I2 + J~ 

+ 2 K z+ s lzi ~S 2 + + 2JC z _~s lzi ~s 2 - + 2 IC +z si + S 2 z 2 + 2 K_ z si_S 2 Z2 

+ (J 7 ++ + fc++) S1+S2+ + (1T4— + /C H —) Si + S2_ 

+ (i7”-++/C.—f) si_S2+ + (j7”—+ /C—) si—S 2 —• (19) 


Here, the parameters and /C are combinations of exchange parameters Jkqk'q' from 
Eq. (J6]) and products of Clebsch-Gordan coefficients: 


J11 


Ji± 


J±i 

J±± 


J±t 


IC ZZ 

1Cz± 

1C ±Z 

1C±± 


/C ±T 


even /~<J 2 m 2 

E rr '~ y JimikO^ J 2 m 2 k'0 

JkOk 

kk’ 


L 'JlJlfcO U J 2 J2fc'0 


(-iJimi -'-(J 2 m 2 

T im Ir-Ci —' _m -• 


1 4- ( — 1 1 2m 2 ^ ven 

rr* _ x ' V \ ' rr JimikO J 2 —m2k'2m2 

Jl^f — ^ / J Jk0k'±2m 2 ~ 

kk’ 


fiJiJi riJ2J2 
U JlJlfcO L/ J 2 J2fc'0 


1 . / -1 \2mi even ^-yj 2 ?7i 2 

,7-* ^ ~t \ -*■) \ ' ,7- Ji— mik2mi J 2 m 2 fc'0 

>4qi/ 4 / ^ i4fc±2mifc'0" 


kk' 


kk' 

X _)_ ^_||2(mi+m 2 ) 

2 

^_-Q2(mi+m 2 ) 

2 

odd ryJimi -'-yj 2 m 2 

E -7- '~ y J\m\kO'^' J 2 ?n 2 fc'0 

JkOk 

kk' 


^ ^ 7Jk±2mik' ±2m 2 


/^yJ 2 J 2 

Ly JlJlfc0 L 'J 2 J 2 fe'0 


ryj 2 m,2 

^ J\— mik 2 mi J 2 —m 2 k' 2 m 2 
riJiJi r^J 2 J 2 
^ JiJikV^ J2J2WQ 


Jfc±2 


!mifc , =p2m 2 " 


kk' 


cJl m l ryj 2 m 2 

^ Ji-mifc2mi l / 2 -m 2 fc'2m 2 

riJiJi (~<J2J2 
u JiJifc0 c 'J 2 J 2 fc / 0 


riJ2J2 ’ 
L 'JiJifc0 u J 2 J 2 fc'0 


1 - (-l) 2m 2 

“ 35 2 


odd 


^C-r — ±-^ JkOk' 

kk' 


±2777-2 " 


iJ 2 m 2 

^ Jim 1&0^J2—777-22777-2 

J 2 J 2 

^J 1 J 1 k0^J 2 J2k , 0 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 


1 / -1 )2mi odd /~i J 2 m 2 

1^-* - V -V \ ' -7- Ji — mi fc2 miJ 2 ?n 2 fc '0 

— ±--- / ; Jl:±2miW 


kk' 


^yJ 2 J 2 

L 'JiJifc0 , - / J 2 J 2 fc'0 


/c 




odd 


Cf' 


C± 


777 2 


1 + ( — l) 2 ( mi + m a) ^ ^J x -m.]k2mx^J2-m 2 k'2m2 

c* / j *^fe±2777lfc / ±27772 ' 


fcfc' 


r^J2J2 

^J 1 J 1 k0^J 2 J 2 k , 0 


tc 


=F± 


1 1 f -i\2(mi+m 2 ) odd syJimi J 2 m 2 

d “r l \ ' - 7 - 0 Ji—mik2mi J 2 — m 2 fc'277i 2 7^)n^ 

-X- 2^ Jk±2m 1 k' Z f2m 2 -^. 7l J, ^.7,J 0 -■ 


(27) 


(28) 


kk' 


°JiJlfc0 U J 2 J 2 fc , 0 


The above relations between exchange parameters are derived using the time reversal 
symmetry [10] and the symmetry relations for Clebsch-Gordan coefficients [25j (see 
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The effective Hamiltonian (19) contains the constant part (the first term), the zero- 


field splitting (ZFS) part (terms 3-6) and the exchange part (the rest of the terms). With 
above equations we are able to identify the situations when J and /C become zero [^| 
The overview of such situations in function of the parity of electrons number on sites 
(N\ and N 2 ) and of the relation between 2m l {i = 1, 2) and g max , is given in 

Appendix A.2| This information is sufficient to find general conditions under which the 


exchange Hamiltonian for two interacting doublets (19) becomes of Ising type. 


3.2. Interacting axial doublet and isotropic spin 

In the case of a doublet (|d]) interacting with an isotropic spin, the derivations are similar. 
Now, we only need to consider the matrix elements of operators referring to the first 
site in Eq. (J7|): 

H efr = PiHPi. (30) 

The matrix elements of I I (JtJ) are given as 




kqq' 


Q ’ 


and the projected exchange Hamiltonian looks as follows: 

i 

Heff ^ ' pXAq'Sl 2l T T KL—q'S\ — \ 

q '=-1 


So 


(31) 


(32) 


Here, spherical components for S 2 are used pa and the expressions for parameters /C 
are given as 


odd 


_ \ ' rr '“'Jim.ifcO 

l^zq' / y «-/fcOl 


1 r J 1 J 1 ’ 
k 

odd 

^ Ji— m\k2mi 


AAq' i ^ ^ Uk : 2 


\m\ 1 q'' 


rih J\ 


(33) 


(34) 


As in the previous case, Eq. (|19J), the exchange Hamiltonian (32) can become of Ising 
type under certain conditions. 


4. Two types of Ising exchange interaction 

With the use of the properties of J and /C in Sec. [3j we discuss the form of the effective 
pseudospin Hamiltonian. Based on the analysis, we show that there are two types of 
the Ising Hamiltonian. 

+ The ZFS terms only appear when at least one metal ion is of non-Kramers type (integer J\ and J 2 ) 
and look artifactual (are not invariant under time reversal) because we describe non-Kramers ions by 
pseudospins s = 1/2. This drawback disappear when the description of such ions is done with integer 
pseudospins, e.g., s = 1 [5B’. This is not necessary here because the exchange part (which is only of 
interest here) is described correctly already within such a simplified treatment. 
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Consider as example a situation when both N± and are odd. For m\ > q m ax /2 


and m 2 > q' max / 2, we obtain all parameters in Eq. (19) zero except JC ZZ , resulting in 

Hj 4/C zz Si Zl S 2 Z 2 i (35) 

i.e. in an Ising Hamiltonian of coaxial type ([Tj) , with an ordering of magnetic moments 
shown in the left plot of Fig. [jja). 

If we diminish the axiality of site 2 so that m\ > q 
reduces to 


c/2 and ra 2 < q' ma J 2, Eq. (19) 


Heff ^^zz^lzi^2z2 C ^ l ^ , z-\-^lzi^2-\- 2/C^— S\ Zl S 2 — 5 


(36) 


where all terms contain S\ Zx . The expression multiplying the latter is a combination of 
pseudospin operators of the second site which can be written in the form 


_ _ &(/C z+ ). 

s 2z' s 2x 2 


9?(/C. 


Z-\- / 


1C' 


S 2y 2 + 




1C' 


■®2z2 j 


1C' = \ IQ 


+ \1C Z+ \ Z . (37) 

Since the coefficients in the first equation are real and normalized to unity, the 
corresponding combination of pseudospin projections can be viewed as a rotated 
pseudospin from initial direction 2 2 towards a new direction z’. Accordingly, the 
coefficients are directional cosines of z’ axis in the coordinate system of site 2. Despite 
the fact that pseudospin is not related to any physical angular momentum but only 
defined via Pauli matrices in the basis of two states Q, one can still define its rotation 
via induced transformations of doublet wave functions ED- Then the eigenfunctions of 
the new pseudospin operator s 2 y will be the same functions Q defined with respect to 


rotated quantization axis z', i.e., |Jj ± m,) / . Then Eq. (36) is rewritten as 

Hu = 41C' zz ,Si Z[ S'2 Z ', (38) 

which has the form of an Ising Hamiltonian of non-coaxial type because the 
magnetization at the second site will not be directed along the corresponding main 
magnetic axis Z 2 but along a different axis z' (left plot in Fig. [ljb)). This makes it 
qualitatively different from the coaxial Ising interaction described by Eqs. 0, f§. 
Accordingly, the two Ising interactions will be called hereafter of type I and II. 

On the other hand, the interaction of a doublet Q with an isotropic spin, Eq. 


(32), can only become of Ising type if both parameters JC +q > and /C_ g / are zero. This is 
achieved when rri\ > q max /2 , in which case the Hamiltonian (32) becomes of the form: 

S 2 *' 


H„ =2 IQ. 


'S l 


zi 


S2 


1C ' ZZ / — \Z\lC z o\ 2 + |/Cz+i| 2 + |/C 2 _i| 2 , 


(39) 


in which the combination 

l 

s^= E 
<?'=-! 


lC zq i S 2g / 
1C', 


(40) 
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Figure 2. (Color online) Examples of symmetric superexchange interaction: (a) via 
p a orbital of a single bridging atom, connecting only m = 0 orbitals on lanthanide 
sites; (b) via the highest occupied orbital of the Nj - bridge j2S], promoting interaction 
between m = ±1 orbitals on lanthanide sites; (c) unpaired electron (S = 1/2) at N| _ 
|28| interacting only with m = ±2 orbitals of lanthanide due to symmetry restrictions. 


is a rotated spin S 2 2 from direction Z\ to z'. This means that the Hamiltonian (39) is of 
non-coaxial Ising type, in which the magnetic moment on site 2 points into a direction 
different from the main magnetic axis on site 1, as indicated in the right plot of Fig. 
[Tfb) (cf. Fig. [TJ a)). Note that the Ising interaction of type I is only achieved when JC z ±i 
in Eq. (40) are both zero, which requires special symmetry of the exchange bridge. 


Table 1. Exchange interaction between two doublets and between a doublet and an 
isotropic spin. I and II denote the two types of Ising exchange interaction, Eqs. (JT|) , 
(35) and Eqs. (381, (39), respectively; N stands for non-Ising exchange interaction. 


Ay 


odd 

mi > <?max/2 mi 

— (?max/2 

even 

mi > <3 Wx/ 2 mi 

— *?max/2 

odd 

"2 > <?max/ 2 

I 

ii 

I 

N 

odd 

m 2 < (/max/ 2 

II 

N 

II 

N 

even 

™2 > (/max/ 2 

I 

II 

I 

N 

even 

m 2 < (/max / 2 

N 

N 

N 

N 

isotropic spin 

II 

N 

II 

N 


The classification for all cases are shown in Appendix B| and Table [l] summarizes 
our main result. It shows that rule 2) (Ising interaction of type I) is not satisfied in 
most cases. We can see also that the type of resulting exchange Hamiltonian merely 
depends on the relative values of mi, m- 2 and g max /2, q' max / 2, respectively. Then we can 
generalize the results in Table [I] over arbitrary doublets (J2]) if by m\ and m 2 we will 
understand not particular doublets Q but the minimal absolute values of the index m 
in the expansion of the corresponding doublet wave function | M) in Eq. ([2]). 
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Table [I] shows that rule 1) does not hold either. Instead, the realization of Ising exchange 
interaction depends on the value of g max (g' nax ) which is determined by the smaller of the 
quantities A max + 1 and fc max (Eq. ([9])). A max is small only for sufficiently high symmetry 
of the exchange bridge. For instance, in the case of linear bridge (Fig. |2](a)) it can have 
the smallest possible value, A max = 0. Its value increases fast with lowering of the 
symmetry of the bridge, and the accompanying increase of the mixing of m orbitals, 
being 2 for the case shown in Fig. [2](b) and 4 for the case shown in Fig. [2|c). When 
A mL on sites become sufficiently large, g max (q' max ) is determined by k max (k' max ), Eq. 
(|8j) , i.e., the smaller of the quantities 2/* + l (= 7 for / orbitals) and 2 J* on corresponding 
sites. Designing a strongly axial crystal field on sites one can achieve m l values as high 
as Ji pj. Then such ions with J,; > 7/2 will a priori provide Ising exchange interaction 
with any other magnetic site according to Table [l] Contrary to what is stated by the 
rule 1), these are precisely the ones which can be called Ising ions. They are listed in 
Table [2] 


Table 2. Ising metal ions. Only valencies of metal ions in known complexes are listed 

mm- 


. f N 

M n+ 

/ 2 

/ 3 

/ 4 

/ 8 

f 9 

f w 

f 11 f 12 

Ln 3+ 

p r 3 + 

Nd 3+ 

Pm 3+ 

X b 3 + 

Dy 3 + 

Ho 3+ 

Er 3 + Tm 3+ 

Ac 2+ 







Es 2+ 

Ac 3+ 


U 3 + 

Np 3 + 

Bk 3 + 

Cf 3+ 

Es 3+ 


Ac 4+ 

U 4 + 

Np 4 + 

Pu 4+ 

Cf 4 + 




Ac 5+ 

Np 5+ 

Pu s+ 






Ac 6+ 

p u 6+ 








6. Conclusions 

In this work, we fold necessary and sufficient conditions to achieve Ising exchange 
interaction in materials involving lanthanide and actinide metal ions. Unless the 
symmetry of the exchange bridge is not very high, these conditions are determined solely 
by the electronic properties of individual metal ions. According to these conditions, by 
far not all metal ions can display Ising exchange interaction (Table [2|. We also fold 
that two types of Ising exchange interaction can arise in these systems, the coaxial 
and non-coaxial ones. In particular, it is established that the Ising interaction between 
an anisotropic doublet and an isotropic spin is of second type. The basic property of 
Ising coupling is the lack of dynamics in resulting exchange doublets, also in applied 
magnetic field. This condition is indispensable for achieving single-molecule magnets 
since it ensures the quenching of quantum tunneling of magnetization in individual 
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exchange doublets. The insight gained in this work will contribute to purposeful design 
of lanthanide and actinide based materials. 

Note added. The Ising exchange interaction between strongly axial ions (Dy, Ho) 
has been recently addressed in the preprint {30] within a different approach from the 
present work. 


Appendix A. Exchange parameters for the effective model 


The relations between exchange parameters, (20)-(29), are derived here. In the 


derivation, we use the property of the exchange parameter under time inversion 

v = (-1 ) q+q 'Jk-,^, (A.l) 

and the property of the Clebsch-Gordan coefficients ra. 


r< c 7 — ( 1 \ b+/3 . / 1 r^a-a 

° aab/3 ~ I--U V 2a + 


(A.2) 


Appendix A.l. Exchange parameters of H e fr 


The relations between the exchange parameters of H e ff, J and /C (20)-(29), are derived 


We obtain the form of J and /C expanding Eq. (18). Among the exchange parameters, 


Jn (20) and K zz (25) are obtained by this direct calculation, whereas for the others, 


additional calculations with Eqs. (A.l) and (A.2) are required. 


Ji± and J ±/, Eqs. (21), (22): 


Using Eq. (A.l), 


Jl± — J fcOfe' 


±2m 2 " 


kk' 


^J 2 m 2 

^ J\ m i kO ^ J2—m2 k' 2 m 2 
ryJiJi 


even ✓'yj 2 ra 2 

\d=2m 2 \ A ( SJ ^ J 2 —m 2 /c / 2m 2 




kk’ 


/- 1 J 2 J 2 

'~ / JlJlfc0'~'J 2 J 2 fc'0 


= (-1 y m2 jj* 


(A.3) 

Since m 2 is integer or half-integer, (—l) 2 ” 12 = (—l)^ 2rri2 . On the other hand, 


applying Eq. (A.2) to site 2, 


Jl± = J k0k , ± 2m 2 

kk' 

= (-1) 2 " 12 Jl±. 

Since k! is even, (—l) fc/ = 1. Therefore, 

Ji± = (-1 ) 2m2 Ji± = 


r'fJimi ^-yj 2 m 2 

^ Ji mi kO ^ J 2 —m 2 k' 2m 2 
riJiJi 

{ ^J 1 J 1 k0^J 2 J2k'0 


(A.4) 

(A.5) 


and we obtain Eq. (21). Similarly, we obtain Eq. (22). 
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J±±, Eq. (23): 


Using Eq. (A.l), 


J±± — E 'Jk± c lm\k' 


'±2m 2 ' 


r^J 2 m 2 

^ J\— mik2mi ^ J 2 —m 2 k'2m 2 


kk' 


^J 1 J 1 k0^J 2 J2k , 0 


= (- 1 )^+^) (Jh 


/c=p2mifc / =F2m2 ) 


kk' 


r'iJimi syj 2 m 2 

^^ J\—m\k2m\ J 2 —m 2 k'2m 2 

rtJ\ Ji rtJ 2 J 2 
Ly JiJi/c0°J 2 J2fc / 0 


_ j'_||2(mi+m2) rj 


=F=F' 


(A.6) 


Using Eq. (A.2) for both sites, 

even 

< 7 - _ \ V i \k+2m 1 +k l +2m. 2 rr 

J±± ~ U Jk±2m 1 k'±2m 2 - 

kk' 

— (—i) 2 ( mi+m2 )j7 ±± . 

From these equations, 

J±± = (-1 )2(m 1+ m 2 ) J- ±± = (_l)2(mi+m a )^ j 

and we obtain relation (23). 


/,^J 2 m 2 
^ J\—m\k2m\ J 2 —m 2 fc'2ra 2 


L 'JlJlfcO C 'J 2 J2fc'0 


(A.7) 

(A.8) 


J± T , Eq. (24): 


Using Eq. (A.l), 


^±=F - X] yJk±2m\k‘ 


'=p2m 2 ' 


yoJimi sij 2 m 2 

^ J\— mik2mi ^ J 2 —m 2 k'2m 2 


kk' 


rtJ 2 J 2 

Ji J\ kO J 2 J 2 k'0 


kk' 

_ ^_|j2mi +2m2 




fc=p2mi/c / d=2m2 ) 


/iJimi 

\* ^J\—m\k2m\ <y J 2 —m 2 k'2m 2 


c 


J 2 m 2 


(- 1 J 2 J 2 

Ly JiJifc0 Ly J 2 J 2 fc , 0 


=F±- 


(A.9) 


Using Eq. (A.2) for both sites, 

even 

<7 _ \ V i\fc+ 2 mi+fc'+ 2 m 2 /r 

c/±=f — / J-J Jk±2m 1 k’^2m2 " 

kk’ 

_ /_-|_^2(mi+m 2 )j^ 

Combinig these relations, 

J± T = (_l)2(mi+m 2 ) i7±T = (_i)2(mi+m a ) t7; . ±j 


/-iJimi /-iJ 2 m 2 

^ J\—m\k2m\ J J 2 —m 2 k'2m 2 


rtJ\ Ji s~t<J 2 J 2 
°JiJifc0 Ly J 2 J2/c / 0 


(A. 10) 

(A. 11) 


and we obtain Eq. (24). 


1C Z± , /C ±2 , Eqs. (J26J), (J27J): 
Using Eq. (j A. 1 h , 


odd 


fczi: i ^ ^ jTfcOfc' 


±2m 2 " 


kk' 


A-fJimi /~iJ2iri2 

^ JimikO^ J 2 — m 2 fc'2m 2 

riJiJi 1 J 2 J 2 

° JiJikO^ J 2 J 2 fe'0 


odd 


^~<J 2 m 2 

1 ( i\2m 2 \ ' / \* JimifcO'- - J 2 -m 2 fc'2m 2 

— ± ( — J-J / fUfc0fc / =F2m 2 J 

kk' 

= - (-1 ) 2ma /c: T . 


flJlJl S~lJ 2 J 2 
^JiJik0^J 2 J 2 k'0 


(A.12) 
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Using Eq. ( |A.2| ) for site 2, 

odd Jimi 

IT — i \ V i\k'+ 2 ni 2 st L 'JimifccA'J 2 -m 2 fc'2m 2 

^z± ~ ± 2_^\ -*■) Jk0k'±2m2 r< J 1 J 1 n J 2 J2 

kk 1 L 'JlJlfc0 L/ J 2 J2fc'0 

= - (-1 ) 2ma /Q ± . 

Here, we used that k' is odd. Therefore, we obtain 


ic z± = -(-l y m 2 K z± = -(-l ) 2m3 /c: T , 


(A. 13) 
(A. 14) 


and Eq. (26). Similarly, we obtain Eq, (27). 


1C±±, Eq. (28): 


Using Eq. (A.l), 


odd 


Y' _ \ ^ rr ^Ji—mik2mi K ~ y J2—rri2k'2m2 

^zbdz / J k±2mik'±2m2 ~ 

kk' 


{ ^J 1 J 1 kO^J 2 J2k'0 


odd 


= (-l) 2 (™l+ m 2) ^2 ( J k 


,C J d 


fc=p2mi k'^f2rri2} 


\* ^J\—m\k2m\ K ~ J J 2 —^ 2 ^2m2 


m 2 


kk' 


r^JiJi r^J2J2 
^J 1 J 1 k0^J 2 J2k'0 


_ ^_^2(mi+m2)]£* 


(A. 15) 


Using Eq. (A.2) for both sites, 

odd 


yoJ2?ri2 

\r _ \ ^ ( i \ H2mi+fc , +2m 2 ^ ^ J\—m\k2m\ J 2 —rri2k'2m2 

^±± — / A U 74±2mifc'±2m 2 ^-yJ 2 J 2 




«/i &0 w J 2 J 2 k'O 


Therefore, 


= (-l) 2(mi+m2) /C ±± . 


/C±± = (-l) 2(mi+m2) /c ±± = (_i)2( m i+ m 2)/c 


/C ±T , Eq. (29): 


Using Eq. (A.l), 


odd 


— _ ^ 74±2 


kk ' 


mi/c'=p2m2 " 


odd 


TT* 


/-iJ2rri2 

^ J\—m\k2m\ y J 2 —fyi2k'2m2 


(A. 16) 
(A. 17) 


r^JiJ 1 ^^ 2^2 

O JiJi/c0 Ly J 2 J 2 /c / 0 


.cf 2 


Jimi 

_ / 1 \2(mi+m 2 ) \ A / rr \* ^Ji—mik2mi'^J 2 —rri 2 k' 2 m 2 

V U WfcT2mifc'±2m 2 J ^-yJaJa 

fcfc' L 'JiJifc0 ( -'J 2 J 2 fc'0 

= (_i)M™i+™ 2 )r ±i (A.18) 


m 2 


Using Eq. (A.2) for both sites, 

odd 

/C ±T = - J](-l) fc+2TO1+fe,+2m2 Jfc± 2 . 


!miA; / =p2m2" 


/^iJ2'm2 

^ Ji—mik2mi^ J2—m2k'2m2 


kk' 


J 2 J 2 

Ly JiJifc0 L/ J 2 J 2 fc , 0 


Therefore, 


= (-l) 2mi+2m2 /C ±T . 


/c ±T = (-i) 2 (mi+m 2 ) /c ±T = (-i) 2 (mi+m 2 ) /c; ± , 


(A.19) 
(A.20) 


and we obtain Eq. (29). 
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Some of J and /C become zero due to (i) time reversal symmetry and (ii) the structure 
of the exchange interaction. 


(i) Time reversal symmetry: 

When magnetic site has odd (even) number of electrons, the total angular 
momentum and its projection m, are half-integer (integer), and thus, (—l) 2mi = —1 
(= 1). Due to the change in sign, some of J and /C (20)-(29) become zero. The 
results are summarized as follows: 


• Jj± = 0 when N 2 is odd, 

• J±i = 0 when N x is odd, 

• K, z ± = 0 when N 2 is even, 

• K-±z — 0 when TV) is even, 

• J±± = J± T = 0 when one of N\ and N 2 is odd, 

• K,±± = /C ±T = 0 when one of Ni and N 2 is even, 

where N { and N 2 are the numbers of electrons on sites 1 and 2, respectively. 

(ii) Structure of the exchange Hamiltonian: 

The values of the exchange parameters J and K, depend on g max 0- For large 
rrii > gmax/2, some of J and /C are zero because either Clebsch-Gordan coefficient 
Cj-mk 2 m or Jkqk'q' becomes zero. 


First, we consider non-symmetric system where g max corresponds to k max . In this 
case, J and /C become zero only when = 0, which holds when k > 2m. 

Since maximum of k is fixed, we have to further consider two cases regarding J: 

(a) dj <: " 1/ 1/2 and k max (/max 2<7j. 

Since 2m < k max is always satisfied, Clebsch-Gordan coefficients are in general 
nonzero, Cj™ mk 2m 7^ 0- Therefore, all J and /C can be nonzero due to the 
reason different from time reversal symmetry. 

(b) Ji > k + 1/2 and fc max = g max = 2{l t + 1/2): 

It is possible to have m which is larger than /c max /2. In this case, we obtain 
Cj m .mk uur , 2 rn = O’ an d J an d /C which include the Clebsch-Gordan coefficient 
become zero. 

The case of N = 6 (J = 0) is not considered since the ion is nonmagnetic. The 
above discussion is summarized as follows: 

• J ±I = J±± = J± T = K± z = 1C±± = JC± T = 0 when mi > h + 1/2, 

• Ji± = J±± = — fcz± = /C±± = ~ 0 w h en m 2 > l 2 + 1/2. 

In symmetric systems, g max < /c max , and the condition changes. Contrary to the 
non-symmetric case, J and /C become zero for m such that fc max > 2m > g max > 0, 

because Cj™ mk2m is not zero for k > 2m, while Jk 2 mk’q' = 0. The selection rule 

becomes 

• J±i = J±± = J± T = JC ±Z = /C ±± = JC ±Zf = 0 when mi > g max /2, 
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• Jr± = J±± = J± T = fcz± = /C±± = = 0 when m 2 > q' ma j2. 

Appendix B. Classification of the effective Hamiltonians 

With the use of the conditions given in Sec. |Appendix X.2\ we study the structure 


of the effective Hamiltonian, (19) and (32), for all cases. The effective Hamiltonian is 


classified into three types, Ising I, Ising II, and non-Ising. Ising I does not change the 
magnetic axes due to the exchange interaction, while Ising II does. 


Appendix B.l. Doublets with unquenched orbital momentum 


The form of the effective Hamiltonian between anisotropic doublets (19) becomes as 


follows (the condition for non-symmetric case is obtained by replacing g max by l + 1/2): 

(i) Both of N\ and N 2 are odd (mi and m 2 are half-integer). 

(a) mi > g max /2 and m 2 > q' ma j2 (Ising, I) 

$lzi S 2 z 2 


H e fi — Jnhh + A 2 

si s 2 

(b) mi > g max /2 and m 2 < q' ma , x /2 (Ising, II) 


(B.l) 


H eff = Jnhh + + K z+ S -^~s 2+ + IC z J^~s 2 - (B.2) 


— Jnhh + IC' ZZ ' 


Si s 2 
Si 21 s 2z ' 


s 1 


Si 


Si s 2 


where 


K„, = At + I^+I 2 . 


. _X(K,+). 

S 2z ' ^ S 2x2 


S(/C 


2+y 


K' 

As As ZZ 

(c) mi < g max /2 and m 2 > q' ma J2 (Ising, II) 
H e ff = Jni\h + fczz 

= J 11 IJ 2 + Ab 


S 2 y 2 + 




K,' 


s 2z2 


Si s 2 
Slz> S 2z2 
si s 2 


+ A +2 Si + — zr^- + K,- Z S\- 


S 2 


S 2Z2 

S 2 


where 




S\z' ~ 


M(1C +Z ) 


$lxi 


%(1C +Z ) 


1 'S '1 1X 1 \ri 

A^ z > z A^ z > z 

(d) mi < g max /2 and m 2 < g^/2 (non-Ising) 

^12:1 ^ 2 z 2 


Slyi + 


/c 2 


K' Z 




H e ff — JllI\I 2 + 

S\z\ 


Si s 2 


Xv | _ ' S, + K.~- z ' A, 

Si Si 

+ /C +J Si + ^ + K_«Si_ S211 


S2 


S2 


(B-3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B.8) 

(B.9) 
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+ K ++ S\ + S2+ + JC-\ —S 1 +S 2 - 
+ /C—|-Sl—S 2 + + K. -Si_S2- 

i'T i r I Y' ^2^2 

— Jn 1 i 1 2 + r^zz— -— 


(B.10) 


-Si S 2 


+ K(/C z+ )^^ - 

S1S2 S1S2 

+ - s(k: +2 )^4^ 

S1S2 S1S2 

5J(/C_|__|_) S\ xl S2x2 ^1j/i^2j/2 


s lzi s 2j/ 2 


2 

sis 2 

3 (/c ++ ; 

^ ^1x1^21/2 ~h S\y 1 S2x2 

2 

-S 1 -S 2 

9 f?(/c+_; 

^ S\x\S 2X2 S\y 1 S2y2 

2 

S 1 S 2 

3(/C + _) 

^ Slxi^2y2 ^l2/i^2rr2 


2 SlS 2 

(ii) TV, is even and N 2 is odd (mi is integer, m 2 is half-integer). 

(a) mi > g max /2 and m 2 > (/ r ' Tiax /2 (Ising, I) 

Heir = Jiihh + IC ZZ S -^ S -^. 

-Sl -S2 

(b) mi > g max /2 and m 2 < q’ mSi J 2 (Ising, II) 

H ef r = J///1/2 + + /C 2+ ^s 2+ + /C 2 _^s 2 _ 

-Sl -S2 -Sl -Si 


(B.ll) 


_ st T T 1 Y' t ^lz\ $2z' 

— J II+l-L‘2 1 K, - , 


Si s 2 


where 


Kz’ = sJ£ 2 zz+\fcz+\\ 


~ _ 3 ft(/C*+) _ 

®2z' ^212 

/V/ / 


9(/C* + ) s , /C 22 s 

j^v s 2j/ 2 “I" t-/ -S2Z2- 

?■ y' 7 7' 


(c) mi < g max /2 and m 2 > g( nax /2 (non-Ising) 


Hoff — Jill 1/2 + + J + jS\ + l2 + J-lS\-l 2 - 


-Sl -S2 

(d) mi < g max /2 and m 2 < g(n ax /2 (non-Ising) 

H cff = Jiihh + IC ZZ S -^ S -^ 

Sl S 2 


(B.12) 

(B. 13) 
(B. 14) 

(B-15) 
(B.16) 

(B. 17) 


J+IS 1 +I 2 + J^lSi-l2lC z+ —J-S2+ + /C 2 _3— !_s 2 _. 

Si Si 


(B.18) 


(iii) Both of N\ and N 2 are even (mi and m2 are integer), 
(a) mi > g max /2 and m 2 > q' max /2 (Ising, I) 

Heff ^Jnhh + icj-^ 8 ^. 

-Sl -S2 


(B.19) 
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(b) mi > g max /2 and m 2 < q' max /2 (non-Ising) 

H eff = JnlJ -2 + + J I+ h~s 2+ + Ji-hh- (B.20) 

Si s 2 

(c) mi < q max /2 and m 2 > q' max /2 (non-Ising) 

H eff = Jnhh + + J+ih+h + J-rh-h. (B.21) 

Si s 2 

(d) mi < q max /2 and m 2 < q' max /2 (non-Ising) 

Heff = Jnhh + 

Si s 2 

+ Ji+Iis 2 + + Ji-hs 2 - + J7+/Si + / 2 + J71t-si_/ 2 

+ J ++ Si + S 2+ + Ji —Si_|_S 2 _ 

+ J7d+Si_s 2 4_ + J — Si_s 2 _. (B.22) 


Appendix B.2. Doublet and isotropic spin 


The effective Hamiltonian between the anisotropic doublet and isotropic spin (32) 
reduces as follows: 


(i) mi > q m ax/2 (Ising, II) 


where 


H e ff ^ ^ !C Z q J 


SlZ! S 


2 q' 


q '=-1 


Si S 2 


y~r Si Zl S 2z , 

— ^zz'~ ^T“ > 

Si >S 2 


1C ', = 


x £ I K 


2 

zq' I ) 


lC zq ' S 2g / 


S *>' = E K- 

g' = -l 2^' 


(ii) mi < g max /2 (non-Ising) 


i 

H 0 ff = 

q '=-1 L 


lC Z q'—~ + ^+q'Sl+ + IC—qtSi— 
Si 


£>2q' 

~S-2 


(B.23) 

(B.24) 

(B.25) 

(B.26) 

(B.27) 
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